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1. BACKGROUND 

In order to understand the distribution and evolution of marine stratocumulus, we 
need to determine how variables such as radiation divergence, surface heat and moisture 
fluxes, wind shear, subsidence, and lapse rate control the dynamical and microphysical 
structure of the clouds through processes such as entrainment and mixing, and in turn, 
how the cloud structure affects these variables. A key to studying entrainment is to be able 
to resolve the structure of the cloud at and near cloud-top. Ideally, we would like to be able 
to resolve both the motion field as well as some tracer that allows us to distinguish well- 
mixed boundary-layer air from that originating in the free atmosphere and entrained into 
the boundary layer. At the top of the stratiform-capped boundary layer, the cloud itself 
can be a tracer of boundary-layer air. Entrained parcels, or regions of mixed boundary- 
layer and free air are clear. The NCAR Airborne Infrared Lidar (NAILS, Schwiesow, 1987) 
is a tool for distinguishing the interface between the clear and cloudy parcels. 

Some progress has already been made in studying the characteristics of the en- 
trainment region of marine stratocumulus. Jensen and Lenschow (1978), for example, 
investigated the sizes of turbulent eddies at the top of a cloud-capped mixed layer using 
data from an aircraft flight leg at cloud-top. They found two peaks in the vertical ve- 
locity spectrum at cloud- top — one with a wavelength of about 1.5 times the depth of the 
mixed layer, and the other with a wavelength of 200 to 300 m that likely was associated 
with entrainment and mixing near cloud-top. Mahrt and Paumier (1982) looked in more 
detail at the entrainment process at cloud top for the same experiment, and estimated a 
scale of 100 to 200 m for the eddies involved in cloud-top entrainment. These results have 
been confirmed in subsequent observational studies by Caughey and Kitchen (1984) for 
nocturnal stratocumulus over land using tethered balloon observations. 

A theoretical analysis of the preferred wavelengths and frequencies of vertical veloc- 
ity fluctuations in the inversion layer capping a clear-air mixed layer and the interactions 
of waves and turbulence at the top of the mixed layer has been carried out by Carruthers 
and Hunt (1985 and 1986), and related numerical simulations are reported by Carruthers 
and Moeng (1987). Again, these analyses suggest vertical velocity wavelengths of a few 
hundred meters in the vicinity of the mixed-layer top. From the standpoint that clouds 
offer a visible tracer of the interface between the free atmosphere and the boundary layer, 
the cloud-capped boundary layer has an advantage over a clear boundary layer. [We note, 
however, that occasionally a thin cloud-free layer has been observed just below the capping 
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inversion (e.g. Rogers and Telford, 1986). This situation, although possibly not common, 
warrants further study to see if differences in cloud-top structure are also detectable with 
NAILS; we recall, on the basis of our inflight observations, that one or more cases of this 
situation exist in the FIRE data set.] 

Thus far, most of the observational evidence for studies of the eddy size and struc- 
ture in the vicinity of cloud top has come from direct measurements of vertical velocity, 
temperature and cloud droplet size distributions. Recently, however, Boers and Spinhirne 
(1987) have reported lidar observations of cloud top which illustrate some of the poten- 
tial that this technique has to offer. They show, for example, a frequency distribution 
and spectrum of cloud-top height for a marine stratocumulus case, and indicate that in 
situ measurements in cloud were obtained concurrently on a second airplane; the in situ 
observations are discussed in Boers and Betts (1988). 

In addition to investigations in the vicinity of cloud top, the NAILS was also pointed 
up at cloud base while flying below the stratocumulus. In this way, it is possible to look 
at cloud-base structure to see how it correlates with measurements in the lower boundary 
layer. A couple of interesting and important questions that can be addressed are: is the 
scale of cloud organization seen in cloud-base height detectable in variables measured at 
aircraft height? How well can the lidar detect multiple cloud layers which, as discussed 
by Nicholls (1984) and Hanson et al. (1988), are indicative of decoupling of the main 
cloud layer from the boundary layer, and if detected, can they also be related to variables 
measured at aircraft height? As Nicholls (1984) points out, this decoupling can make a sig- 
nificant difference in the subsequent evolution of the cloud cover, and therefore significantly 
affect the radiation budget of the cloud-capped mixed layer. 

2. NAILS 

The NCAR airborne infrared lidar system (NAILS) used in the 1987 FIRE exper- 
iment off the coast of California is a 10.6 /zm wavelength CO 2 lidar system constructed 
by Ron Schwiesow and co-workers at NCAR/RAF. The lidar is particularly well suited for 
detailed observations of cloud shape; i.e. height of cloud top (when flying above cloud and 
looking down) and cloud base (when flying below cloud and looking up) along the flight 
path. Schwiesow et al. (1988) have presented some examples of NAILS lidar data from 
FIRE. 

A detailed discussion of the lidar design characteristics is given by Schwiesow (1987), 
and only a brief summary is listed here (Table 1). The lidar height resolution of ± 3 
m allows for the distance between the aircraft and cloud edge to be determined with 
this accuracy; however, the duration of the emitted pulse is approximately 3 /is, which 
corresponds to a 500 m pulse length. Therefore, variations in backscatter intensities within 
the clouds can normally not be resolved. Hence the main parameter obtainable from the 
lidar is distance to cloud; in some cases we may also be able to determine the cloud depth. 

During FIRE the lidar was operational on 7 of the 10 Electra flights, and data were 
taken when the distance between cloud and aircraft (minimum range) was at least 500 m. 
The lidar was usually operated at 8 Hz, which at a flight speed of 100 ms -1 translates 
into a horizontal resolution of about 12 m. 
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TABLE 1 

NAILS Lidar Characteristics (modified from Schwiesow, 1987) 


Operating wavelength 
Beam width 

Lidar return sampling rate 
Height resolution 
Track resolution (variable*) 
Design backscatter coefficient 
Maximum range from aerosols 
Maximum range from cloud 
Minimum range (dead zone) 
Signal channels (polarization) 


10.6 fim 

16 xlO -5 radians 
100 MHz 
±3m 
~ 10 m 

5 x 10- # (m" 1 sr- 1 ) 
2 km 
5 km 
500 m 
2 


'Depending upon pulse frequency and flight speed. 


The backscatter as function of time (equivalent to distance) for each laser pulse 
is stored in digital form on magnetic tape. Currently, three independent variables are 
available to the investigators on the FIRE Electra data tapes: lidar range to cloud, strength 
of return (relative power) and pulse width of return, which is related to penetration depth. 

The values of these variables should currently be interpreted^with some caution. In 
the future we hope (i) to test the detection algorithm in order to insure that the actual 
cloud edge is detected, (ii) to determine under what circumstances the lidar is able to 
detect the heights of both cloud base and cloud top, and (iii) to correct the data for pitch 
and roll errors. 

3. EXAMPLES OF NAILS DATA 

The following example illustrates only one of many ways the lidar can be used in 
conjunction with other sensors in order to study stratocumulus clouds. Schwiesow et al. 
(1988) present time series of lidar-derived cloud top height and radiation temperature at 
cloud top from the Barnes PRT-5 sensor. Fig. 1 shows the cloud top height, z c , measured 
by the Electra from 211100 to 211400 PDT on June 30, 1987. The lidar shows numerous 
variations in the cloud top height on scales of hundreds of meters. In the following we 
will examine the association between the variation of cloud top height and cloud top 
temperature. 

In Fig. 1 we have selected a segment of the cloud top for calculating a reference 
value of entropy. If we assume a moist adiabatic temperature gradient in the cloud, then we 
can calculate the temperature of the cloud surface based on the lidar measured cloud top 
height. This lidar-derived cloud top temperature, Tj, can then be compared with the PRT- 
5 temperature, T r ; see Fig 1, center panel. Differences between these two, AT = T r — Tj, 
are shown in the bottom panel of Fig. 1 and they can be due to several effects: 
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Fig. 1. Top panel: Lidar-derived cloud top height (z c ) averaged to 1 s -1 . Middle 
panel: PET-5 radiation temperature (T r , bold line) and lidar-derived cloud top temper- 
ature (T|, thin line). See text for explanation of arrows. Bottom panel: Temperature 
difference, AT = T r — Tj. The thin line shows AT = 0. The distance between tick 
marks on the abcissa is ^ 1 km. 


(1) Radiative cooling and the subsequent formation of cold pools may lead to a pattern 
in which local depressions in the cloud top surface have predominantly negative 
AT. 

(2) Entrainment and the ensuing evaporative cooling may likewise lead to local depres- 
sions in the cloud top surface with negative AT. 

(3) It is not clear to what extent the lidar-derived cloud edge coincides with the region 
which dominated the PRT-5 signal. If the cloud is very thin, then the PRT-5 
temperature may include contributions from a substantial cloud depth, whereas the 
lidar may detect the actual cloud edge. This would lead to regions of positive AT. 

(4) The approximation of moist adiabatic temperature gradient may not be valid for 
the cloud. 

An examination of Fig. 1 does not result in obvious correlations between z c and AT; 
hence without further investigation we can not conclude that “domes” contain warm air, 
whereas depressions between them contain cooled air. Examples of “thin” cloud appear 
likely in at least four segments, as indicated by the arrows in Fig. 1. 
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